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ABSTRACT 

We present high-resolution images of the strongly bipolar planetary nebula NGC 6537, 
obtained with Hubble Space Telescope and with the infrared adaptive optics sys- 
tem on the Very Large Telescope. The central star is detected for the first time. 
Using the multi-band photometry and constraints from the dynamical age of the neb- 
ula, we derive a temperature in the range 1.5-2.5x 10 5 K, a luminosity ~ 1O 3 L , 
and a core mass M c w 0.7-0.9 M . The progenitor mass is probably in the range 
Mi = 3-7 M . The extinction map shows a largely symmetric, and compact dust 
structure, which is most likely a shell, located at the neck of the bipolar flow, only 
2-4 arcsec from the star. The dust shell traces a short-lived phase of very high mass 
loss at the end of the AGB. The dynamical age of the shell and bipolar lobes are 
very similar but the morphologies are very different. The data suggests that the mass 
loss during the ejection of the compact shell was largely spherically symmetric, and 
\ the pronounced bipolar ity formed afterwards. The dynamical ages of the bipolar lobes 

and dust shell are similar, which is consistent with suggestions that bipolar structures 
form in a run-away event at the very last stages of the AGB mass loss. The inner edge 
of the dust shell is ionized, and PAH emission is seen just outside the ionized gas. We 
associate the PAH emission with the photo-dissociation region of the molecular shell. 
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1 INTRODUCTION 

At the end of the asymptotic giant branch (AGB) phase, 
stars experience high mass loss rates (10 -7 -10 -6 Mq yr _1 ), 
which may increase by a further ord er of magnitude dur- 
ing helium shell flashes (Habine 1996). The high mass loss 
terminates once the hydrogen envelope of the star is almost 
fully removed. During the subsequent post- AGB phase, the 
temperature of the central star increases, and eventually the 
surrounding ejected material is ionised, and detected as a 
planetary nebula (PN). 



* Based on observations with European Southern Observatory, 
Very Large Telescope with an instrument, ISAAC and NAOS- 
CONICA (the proposal numbers: 65.D-0395A, 72.D-0766A). 
f Based on observations made with the NASA/ESA Hubble 
Space Telescope, obtained from the ESA/ECF Data Archive. 
HST is operated by the Association of Universities for Research 
in Astronomy, Inc., under NASA contract NASA 5-26555. These 
observations are associated with program no. 8345. 



During the so-called superwind, the objects suffer high 
self-extinction from dust which forms in the eject a; the 
absorbed radiation is re-emitted at infrared wavelengths. 
The IRAS photometric observations show that the infrared 
colour temperature decreases from 1000 K in AGB stars 
to a few hundred Kelvin in post-AGB stars an d PNe 
(|van der Veen and Habindll988MZiilstra et "aLll200lh : as the 
star evolves, the circumstellar shell expands and the tem- 
perature of the circumstellar shell decreases. 

Planetary nebulae commonly show non-spherically 
symmetric morphologies, often elliptical or bipolar. These 
may trace the interaction of a fast wind (> 1000 km s -1 ) 
from the central star of the PN, with the sl ow wind (20- 
30kms _1 ) dating from the AGB towoklll982h . The expan- 
sion velocity of the AGB shell is 20-30 kms -1 while the fast 
wind from the PN is more than 1000 km s _1 . This interaction 
of the winds can enhance shapes which are already present 
in the original AGB wind. However, it is still not clear how 
the original shape is created. There is evidence that bipolar 
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PNe tend to have higher mass progenitors (|Kastner et alJ 
1996), from h elium and nitro gen enrichments characteristic 
of type I PN (IPeimbertlll97^ . These have likely progenitor 
masses of > 2.5 Mq. Magnetic fields have been suggested, or 
a bin ary companion may giv e rise to a circumbinary dusty 
disk (|Balick and Frankll2002f) . A war ped disk h as been sug- 
gested to give rise to multipolar PNe ( Icke 2003). Some AGB 
stars show a cold detached dust shell. Infrared images have 
resolves some of these s hells, which appear not to show devi- 
ation s from symmetry (|lzumiura et alJll99d lOlofsson et alJ 
l200(tl . 

Here, we report high spatial resolution images of the PN 
NGC 6537, which shows one of the most pronounced bipolar 
outflows seen among PNe. It is a ty pe-I PN, which is over - 
abundant nitrogen and helium fc.f. IPottasch et aTll2000bh . 
The cen tral star of NG C 6537 has not previously been de- 
tected (|Pottaschll2000ar ) . We use HST images to create an 
extinction map: this traces the dust at far higher resolution 
than can be achieved using far-infrared dust emission. We 
find that the extinction inside this nebula is not uniform, but 
is concentrated in a compact dust shell. The central star is 
detected for the first time, only 0.5 arcsec away from the 
brightest part of the nebula. Using the photometry in the 
optical HST images and the infrared adaptive optics image, 
the luminosity of the central star is constrained. 



2 OBSERVATIONS 

We observed NGC 6537 with the Adaptive Optics (AO) sys- 
tem, NAOS, and the infrared camera CONIC A, on the Very 
Large Telescope (VLT) on the 6th of March 2004 (in UT). 
The exposure time was 28 minutes. We used the S27 cam- 
era with a pixel scale of 27 milliarcsec. The Ks-band filter 
was used. The dichroic for AO was VIS, i.e., the visual cam- 
era was used. The AO reference star was S 300100217909 
(1/=14.8 mag) from the Guide Star Catalogue, which is ap- 
proximately 22 arcsec away from NGC 6537's central star. 
The approximate spatial resolution of the final image is 
0.22 arcsec (RA orientation) and 0.16 arcsec (Dec orienta- 
tion) , estimated from nearby stars; but the spatial resolution 
varies slightly even within the image. The weather was clear. 
The internal calibration lamp was used for flat fielding. The 
jittering technique was used to minimize the effect of hot 
pixels. The sky background was estimated from a medium 
of jittered frames with slightly different positions. 

NGC 6537 was observed with the Hubble Space Tele- 
scope (HST) in the programme 6502 (P.I. B. Balick). In 
total, HST images were obtained with 8 different filters 
(TableQJ. We retrieved the images from data archives at 
the Space Telescope - European Coordinating Facility (ST- 
ECF). The data were taken on the 12th of September 1997 
with the Wide Field Planetary Camera 2 (WFPC2), which 
has four 800x800 pixel detectors. The main region of in- 
terest for this paper, close to the centre, is covered with 
the Wide Field Cameras which have a 0.1 arcsec pixel scale. 
The spatial resolution is limited by the pixel scale, as the 
point-spread function is under-sampled, but slightly better 
than that because multiple images were taken for a single 
band (with the exception of F487N and F547M). We use 
the pipeline-reduced data. The PHOTFLAM parameter is 
used for the flux calibration. Spike noise due to cosmic rays 



Table 1. The log of HST observations. 



Filter 


Wave 
[A] 


BW 

[A] 


CS 




Lines 


F673N 


6732 


30 


yes 


400secxl,500secxl 


[SII] 


F631N 


6306 


13 


yes 


400secx3 


[OI] 


F658N 


6591 


30 


no 


400secx3 


[Nil] 


F656N 


6564 


54 


no 


400secx3,40secxl 


Ha 


F547M 


5483 


205 


yes 


40sec 




F502N 


5013 


48 


no 


300secx2, 400secxl 


[OIII] 


F487N 


4865 


32 


no 


500sec 


H/3 


F469N 


4695 


17 


no 


400sec 


[Hell] 



Wave: wavelength 

BW: band width (for photometric calibration; not equivalent 
with FWHM) 

CS: detection of the central star 

T ex : exposure time per frame (some bands have multiple expo- 
sures) 

Line: major lines inside the filter 



was removed using the exposures with slightly different posi- 
tions if multiple images were obtained within the same filter, 
exceptions being F487N and F547M which have a single ex- 
posure for each. For F487N, we removed the cosmic ray noise 
by taking the medium over nearby regions. 

Our PAH image was obtained with ISAAC on the VLT. 
The data were acquired on 13th of August, 2000 in service 
mode. We used the eclipse and IDL packages for data reduc- 
tion. Images were obtained with two filters: NB_3.28 (PAH 
band) and NB-3.21 (continuum). Total exposure time was 
7min for each band. Optical seeing was 0.56-0.83 and 0.47- 
0.66 and the measured infrared seeing was 0.53-0.61 and 
0.47-0.54 arcsec (FWHM), respectively, which was measured 
from the nearby stars. After the flux calibration with refer- 
ence to a standard star HR 7446, the continuum was sub- 
tracted, and a PAH image was created. 



3 RESULTS : THE CENTRAL REGION OF 
NGC 6537 

The entire field at Ha is shown in Fig.Q It shows an ex- 
tended bipolar outflow, and a dense core with a size of 
^5x5 arcsec 2 at the centre. The bipolar lobes have the ap- 
pearance of an hourglass, seen close to edge-on; the north- 
western lobe is slightly fainter and may be pointing away 
from us. At the faintest levels, the nebula extends beyond 
the WFPC2 field of view. The lobes extend over 1 arcminute, 
and are oriented roughly NE-SW. The core consists of bright 
arcs tracing a shell surrounding an irregular cavity (Fig. |2| • 
The cavity is between 2.5 and 4 arcsec in diameter, with 
the long axis roughly (but not perfectly) aligned with the 
bipolar flow. The arcs have thickness of approximately 1 arc- 
sec. There is a much fainter extended shell with a diameter 
of 6 x 2.5 arcsec, with the major axis oriented EW. Apart 
from the cavity elongation, none of the structures shows an 
obvious alignment with the bipolar axis. 

The Ks-band image is shown in Fig.|3| Several arcs are 
detected, which are the same as those found in the central 
region of the Ha image, as discussed later. The dominant 
source of emission in the Ks-band is likely to be free-free and 
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Figure 1. The entire field of the planetary nebula NGC6537 in 
Hq. North is top and east is left. The x-axis is the RA offset in 
arcsec and y-axis is the Dec offset in arcsec, for all of the images 
in this paper. The unit is erg cm _2 s _1 arcsec - 
log (I U )=-1A, -13.6, -13.3, -13, -12.6, -12.3, 



2 , contours show 
-12. 



Table 2. Measured fluxes of the central star. The error includes 
only the measurement error of aperture photometry. Systematic 
errors, such as contribution of nebula into the aperture or the 
error of calibration star (for Ks band) are not included. 



Filter 


Wavelength [A] 


Magnitude 


Flux [Jy] 


F673N 


6732 


21.74dz0.54 


[11.0 zb 7.2] x 10~ 6 


F631N 


6306 


21.34±0.33 


[14.2 zb 5.1] x 10- 6 


F547M 


5483 


21.56zb0.22 


[8.6 zb 2.0] x 10- 6 


Ks 


2.18X10 4 


18.63zb0.61 


26 x 10- 6 




-2 




-4 - 



-1 4.0 



Figure 2. The central region of Ha image in 

the unit of ergcm - 2 s _1 arcsec -2 , contours show log 

(J„)=-14, -13.6, -13.3, -13, -12.6, -12.3, -12. (same as 
Fig.P 



bound-free continuum emission from the ionized gas. Only 
regions of higher emission measure in the ionised region are 
detected at the Ks-band. The much fainter bipolar wings are 
not detected. 
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Figure 3. Ks-band image (in linear scale). The lines show the 
position of the central star. 



3.1 Detection of the central star 

The star is clearly detected, separated from the extended 
nebula in the high spatial resolution images with HST and 
NACO. The coordinates of the newly discovered central star 
are RA=18 h 05 m 13!03, Dec=-19?50. / 34 / /4 (J2000). 

The Ks magnitude of the central star is 18.6 =b 
0.6 mag, which is derived from aperture photometry with 
a radius of 5 pixels. The magnitude is referenced to 
nearby stars in the same frame: 2MASS 18051346-1950277 
(Ks=13.382±0.166 mag) and 2MASS 18051378-1950407 
(Ks=12.674±0.040 mag). The uncertainty of 0.61 mag is the 
measurement error only: the uncertainty from the magni- 
tude calibration is not included. There might be some flux 
contribution from the nebula within the aperture, which 
may not have been removed completely. The background 
level (including the contribution of nebula) is estimated from 
circular with inner radius of 150 % and outer radius of 200 
% of the apeature size. 

The central star is also detected in three HST filters, 
two of which are shown in Figs. 21 and |5] The results of aper- 
ture photometry with a radius of 2 pixel are listed in Ta- 
ble^ The largest uncertainty of the magnitudes is due to 
the contribution of nebula and the flux missed outside of 
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Figure 4. HST F673 band image (in logarithmic scale). The lines 
show the position of the central star. 




-4 -2 2 4 

Figure 5. HST F547M band image (in logarithmic scale). The 
lines show the position of the central star. 



the radius. These systematic errors are not counted in these 
errors. F547M is the only medium-band filter and this fil- 
ter is sensitive to continuum. In addition, the central star 
is detected in two narrow band niters. The remaining five 
niters which were used for the observations did not give de- 
tections. The star is too faint to be detected in the narrow 
blue niters (F487N, F469N and F502N; Fig.0, due to ex- 
tinction, while the F656N and F658N niters contain strong 
ionised lines within the filter transmission (Ha and [Nil]), 
and the contamination of these ionized lines from the nebula 
prevents detection of the star (Fig.0. All three HST filters 
where the central star was detected have only weak atomic 
line s within their w avelength coverage. 

IPottaschl fcOOOaT) reported a non-detection of the central 
star, using the same HST data. The under-sampling of the 
Wide Field Camera can make it difficult to distinguish a 
faint star from spike noises due to the cosmic ray hits, and 
this may be the reason for their non-detection. However, 
comparison with the well-sampled Ks-band images confirms 
the central star detect i on in three images. 

Kale r and Jac obv ( 1989) predicted that the ce ntral star 
of NGC 6537 should be about V = 19.67 mag, and IPottaschl 
(2000a|) estimate V = 22.4 mag. We obtained magnitudes 
in a 'continuum band' at 5483 A, which is very close to 
the V-band. Assuming that the visual magnitude is the 
same as that in the F547M band, 21.6±0.22mag and ig- 
noring the systematic errors, our measured magnitu de is 
just between the predict ed values of IPottaschl (|2000a|) and 
(Kaler and Jacobj fiflSS) • 

A dete c tion o f the central star has been suggested by 
iReav et all Jl984), as a 'stellar- like condensation to the 
east of centre.' H owever, iGathier and Pottaschl ([l988) and 
IPottaschl {2000a ) show that this detection is not a star but 
a condensation in the nebula. We agree with the latter view, 
since the nearby condensed part of nebula is only ^0.7 arc- 
sec away from the star, and co uld not have bee n resolved 
from the star at 2 arcsec seeing (Re av et"aLl[l9 84 > ). 

3.2 The extinction map 

The extinction map is estimated from HST Ha (F656N) 
and H/3 (F487N) images. We assume an electron tempera- 
ture of 15 000 K, and a n electron density of 1 x 10 4 cm -3 
(|Pottasch et alJl2000b|). The Ha and H(3 intens ity ratio of 
2.79 (Case B) is used (|Storev and Hummerll995l) . The sensi- 
tivity is limited by 11(3 data (Fig.0. The image is smoothed 
to 0.2 arcsec pixel scale. This also removes any difference in 
resolution between Ha and H/3. 

Fig. |7| shows the extinction map of the core region of 
NGC 6537. The extinction within the core is non-uniform, 
showing the presence of internal extinction within the PN. 
At the centre, the extinction E(Rf3 — Ha) is 1.4-1.5 mag, 
while it increases up to 2.2 mag at 4 arcsec away. The thick- 
ness of the high extinction region is about 2 arcsec. This high 
extinction shell shows the presence of a detached dust shell, 
located far inside the extended planetary nebula. The inner 
cavity of the dust shell shows an elongation along the bipolar 
outflow which is towards the north-east (=left-top in the fig- 
ure), and the south-west (=right-bottom). The central star 
is located at the centre of the elongated dust shell. There 
are small 'patch-like' scale structures inner region the shell 
within 2-3 arcsec from the central star. 
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Figure 6. The central region of H(3 (I u ) image Contour show 
log (J„)=-14, -13.6, -13.3, -13. 




2.2 



Figure 7. The extinction map. The contour lines show 1.4, 
1.6, and 1.8 mag in E(H(3 — Ha), and smoothed to 0.14 arcsec 
pixel -1 . The central star is marked as a cross. At the centre, 
the extinction is 1.4—1.5 mag, while the extinction is higher than 
2 mag at ~4 arcsec from the central star. White points inside 
the extinction map (e.g. RA offset ~ —2 and Dec offset of ~ —4) 
are caused by cosmic ray hits on the detector, in the Hf3 image. 



ICuesta et al.l (|l995f) suggested a reduction in extinction 
at the centre, using their optical spectroscopy across the cen- 
tre region. We confirm this, and show the extinction peaks 
in a shell structure. 

The average e xtinction toward s NGC 65 37 is reported 
as E(B -V) =1.23 (IPottaschll2000aD . A^=3A . (ICuesta et al.1 
1995), or (7=1.79-2.16 (lAshlev and Hvlandl Il988j) . These 
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Figure 8. PAH image in Jy arcsec -1 . 



correspond to E(H0 - Ha) =1.50, 1.33, 1.48-1.79, if the 
inter stellar extinction law (jCardelli. Clavton and Mathisl 
119891 ) and A(V)/E(B - V r )=3.1 are assumed. These val- 
ues of extinction were derived from hydrogen recombina- 
tion lines, presumably weighted towards the lines originating 
from high Ha intensity region. In our map, that region has 
E(Hf3 — Ha) = 1.5 mag, which is within previous measure- 
ment uncertainties. The average E(Hf3 — Ha) of the whole 
region (denned as E(H/3 — Ha) higher than 0.2 and smaller 
than 3.0 mag) is 1.25, and the median is 1.20 mag, which is 
below E(Hf3 — Ha) at the high Ha region. 



3.3 Hierarchy of Ha, PAHs and extinction 

Fig.|H|show the PAH image of the central region. The PAH 
emission is concentrated in an elongated shell. Comparing 
with the Ha distribution (Fig.0 shows that the arcs seen 
in the Ha image are also emitting PAHs, but the shapes are 
not fully coincident: the PAH distribution traces the outer 
edge of the optical arcs. Slices through the Ha, PAH and 
extinction images are plotted in Fig.|5| This shows a radial 
hierarchy in the shells and rings: the Ha rings are found in 
the innermost region, and PAHs are found at and outside of 
the Ha region, where the extinction is E(B — V)~1.8 mag, 
slightly higher than the extinction towards the central star. 
The extinction increases further outside of the PAH emitting 
region. 

To first approximation, the dense Ha emission traces 
the inner edge of the dust shell, while the PAH emission 
could be located in a photo-dissociation region. However, 
confirmation of such a model would require detection and 
mapping of molecular emission. 



4 DISCUSSION 
4.1 The central star 

The magnitudes of the central star allow us to estimate the 
luminosity of the central star, and to discuss its evolutionary 
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Figure 9. A slice through the extinction map, Ha map, and PAH 
image at constant RA (minus is the south) across the central star. 
Extinction is smoothed to 0.8arcsec pixel -1 on image. Bumps at 
Dec offset =—6.2 and 2.5 in the extinction map marked with 
crosses are artificial (one large spike noise spread to other pixels 
because of the smoothing). 



stage. To derive the luminosity, the effective temperature 
and extinction towards the centre are needed. 

For the extinction towards the central star we take the 
central extinction measured from Ha/H/3. This assumes that 
the ionized region is located within the dust shell. The hi- 
erarchy discussed above supports this distribution. We ig- 
nore back- scattering which may make the observed extinc- 
tion slightly higher. If this were important, the dense Ha 
filament located very close to the position of the star could 
be expected to show lower extinction. However, the filament 
shows identical extinction to the general inner region. 

We estimate the Zanstra temperature using our vi- 
sual magnitude of the central star and mean ex tinction, to- 
gethe r with the total flux of H/3 and He 4861 A (|Reav et al.l 
1984). The Zanstra temper ature estimated from helium is 
340,000 K. The estimate of lPottaschl (l2000aTl shows a higher 
Zanstra temperature of 500 000 K. Since we can not estimate 
the systematic error in our measurements of visual magni- 
tude, this higher temperature could not be ruled out totally 
at this point. 

Given the very high temperature, the various uncer- 
tainties on the Zanstra temperatures need to be considered. 
Some cause the Hell Zanstra temperatures to be too low: high 
optic al depth, internal dust extinction, helium overabundance 
iGruenwald and Vie gas 2000). A binary companion could make 
the star appear brighter which would also lead to a lower Zanstra 
temperature. An excess of stellar photons in the extreme ultravi- 
olet, shortward of the He ionizatio n edge, can cause the Za nstra 
temperature to be overestimated. iHenrv and Shipmanl (1986) dis- 
cuss the importance of composition: He-poor atmosphere cause 
the Hell Zanstra temperature to be overestimated, however this 
effect diminishes at the highest temperatures and is perhaps less 
likely for a type-I (He-rich) PN. lGabler et alJ Il99ll) mention that 
shocks within the fast wind can increase the EUV emission and 
lead to overly high Zanstra temperatures. 

From the previous section, the extinction towards the central 
star is E(R(3 - Ha)=1.40, which is equivalent with A(V) = 3.71. 
The magnitudes for the four bands where the star is detected, 
are corrected for this extinction. We fit the resulting fluxes with 
a black body (Tig. 1101 to derive the luminosity of the central star 



(L*) (Table|3J. Initially, parameters in Pottasch (2000a) are used; 
the black body temperature is Tbb=500 000 K, the distance is 
Do=2.4 kpc, and the stellar radius is ro = 1.16 x 10~ 2 Rq. With 
these parameters (case la), the fit is within 1 a in the optical but 
outside of this range at F547M and F631N. 

For the next step, we scale the blackbody by a factor a 2 , i.e. 
the solid angle is Q « 7r(ro/Do x a) 2 . The lowest flux within the 
uncertainty (1 a) is found at a=1.16 (case lb) and highest flux is 
found at a=1.38 (case Ic). The luminosity is calculated with two 
possibilities for each case, adopting the scaling factor for distance 
and for radius, respectively. We assume Tbb is equivalent with 
effective temperature T e ff . Consequently, the luminosity range is 
[7.5-10.1] xlO 3 L (TableEJ. 

Finally, Tbb is varied. A lower estim ate of the te mperature 
of the central star has been reported, and Casas sus et alJ (2000) 
estimated 150 000 K from the ionization levels of the nebula. The 
fitted result corresponds to case (II) in Tabled The luminosity is 
an order of magnitude lower in this case. On the other hand, with 
our estimated Zanstra temperature of 340 000 K, the luminosity 
range is [1.6-4.6] x 10 3 L (Table Case Ilia is for the lower 
limits of the magnitudes errors, and case Illb is for the upper 
limits. The black body lines of case III are not plotted in Fig. llQI 
because the lines overlay with those of case lb and Ic on this 
plotting scale. 

We plot these estimated L* and T e ff (assuming T\ i \ ) =T e ff) 
on the HR diagram with a comparis on with t he theoretical evolu- 
tionary tracks for solar abundances (Blocker 1995). Fig. llll shows 
that the central star is within the highest luminosity as a PN 
and temperature range if T e ^ is 500 000 K, whilst if T e ff is lower 
than that, the central star is already in the white dwarf (WD) 
phase. Different estimates for T e ff give totally different luminosi- 
ties. Nevertheless, the models tend to coincide with the higher 
mass tracks. 

A further constraint on the stellar mass derives from the time 
scales for the post-AGB evolution dGesicki and Ziilstral feoOCft . 
Higher-mass stars evolve much more rapidly to higher tempera- 
ture. All models accelerate when evolving through the knee of the 
HR diagram, and evolution slows dramatically around 10 3 Lq for 
the highest mass tracks, and around 300 Lq for the 3 Mq track. 
Indicative post-AGB ages for different tracks are shown in Fig. 

El 

ICuesta et alJ Jl995l) reported that the central region 
of NGC 6537 has an ex pansion velocity of 18kms _1 . 
ICorradi and Schwarzi {]_993) also found a similar velocity. 
In previous discussion, we obtained distance of 0.9-2.4 kpc. If we 
assume the distance of 1.5 kpc, the dynamical age of the dust 
shell should be about 1600 years, assuming that the shell has 
been constantly expanding with this speed and that the current 
radius is 4arcsec. Models with core masses in the range 0.696 
- 0.943 gives a possible fit for a luminosity in the range 

2.5 < logL < 3. Lower-mass tracks evolve too slowly before 
reaching the dashed region. The best fit track has a core mass of 
~ 0.7 0.9Mq and progenitor mass of 3-7 Mq. 

The comparison of ages and stellar parameters favours a 
lower temperature for the central star of NGC 6537, in the re- 
gion 1.5-2.5xl0 5 K. This is below the Hell Zanstra temperature, 
suggesting the latter may have been overestimated. 

The progenitor mass heavily depends on the mass lost dur- 
ing the AGB phase, and is less well constrained than the core 
mass. The relation between final and initial mass may differ from 
the ones of th e Blocker tracks. However , the derived core mass is 
high for PNe (jGesicki and Zii'l stra 2000) and is consistent with a 
higher-mass pr ogenitor. The n itrogen and helium enrichment of 
this typeT^Nj|Peimber^273) also suggests a high-mass progen- 
itor (c.f. lPottasch et alJl2000b|j . 

The bipolar lobes have velocities of ~ 350kms _1 . This gives 
a similar expansion age to that of the core. This suggest the bipo- 
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Figure 11. T he estimated luminosity and T e ff of the central star of NGC 6537, plotted on the evolutionary tracks of post-AGB stars 
Blocker (1995). The dashed line encloses the estimated temperature and luminosity range (the asterisks show the fit calculations). The 
numbers in the parenthesis are Mzams an d the core mass at the tip of the AGB, in solar units, for iBlockeri (jl995lV s models. The 
numbers on the plus signs show the age in 10 3 years after the AGB phase. The bold lines on the evolutionary tracks indicate the age 
of 1000-2000 years after the AGB phase, suggested by the dynamical age. This age range is estimated from the size of the dust shell. 
M c w 0.7 - 0.9 M and M ZAM s of 3-7 M is likely for NGC 6537. 



lar structure and the dust shell formed near-simultaneously, and 
were present very soon after the onset of the post-AGB evolution. 

4.2 Dust within NGC 6537 

The temperature of the circumstellar dust shell decreases as the 
circumstellar shell expands and the central st ar evolves from AGB 
star to post-AGB star and PN (|van der Veen and Habi^[l988). 



The cool far-infrared excess found in PNe is the remnant of the 
large mass loss phase at the end of AGB phase. 

The temperature of the dust shell is estimated from the cen- 
tral star parameters and the distance, discussed in the previous 
section. The high extinction is mainly found within 2-4 arcsec 
from the star. Ignoring the projection effects, an inner radius of 
2 arcsec and outer radius of 4 arcsec are assumed. We can esti- 
mate the dust temperature from d^fR* = {T^/T^) 2 , if the 
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Figure 10. The extinction-corrected flux of the central star, and 
the black body fits. The parameters are listed in Table|3| Case Ic 
and case II have almost the same line, the dashed line is for case 
Ic, and the blue line is for case II. 



Table 3. Fitting results to the fluxes. 
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Figure 12. The measured infrared to radio flux of the 
core of NGC6537. The data are from our IS AAC ob- 
servations at L and M-ba nds, I S O/SWS, IRAS. |Whitelockl 
fmsl. iPhillms and Cuestal EflSl), |Cohen and Barjfld <mfl, 
and Hoar e et al.l J1992). The solid line is a total of calculated 
free-free and bound-free emission, scaled to the radio data. Two 
black body curves are shown (dashed lines) to emphasize the dust 
excess on top of the extinction corrected free-free and bound-free 
continuum. 



dust grains are directly heated by the central star, where is 
the distance of a dust particle from the central star and R* is 
the stellar radius. For T e ff=3.4xl0 5 K, scale a = 2.6 (Case Ilia 
in Table|3j, is 45 K and 32 K at 2arcsec, and 4arcsec in ra- 
dius, respectively. This can be compared to the dust condensation 
(and formation) temperature, about 1 000 K for silicate. The dust 
grains formed much closer to the star, during the AGB, and the 
grains have moved outwards. 

Infrared spectra (Tig ll21 show two peaks in the dust excess. 
The colder shell has a peak at about 50-60 /im, which corresponds 
to temperature of about 50 K. This is consistent with the above es- 
timated temperature from the dust radius. The hotter dust peaks 
around 10/im, with a temperature around 300 K. This hotter dust 
cannot be located much closer to the star, as the optical and in- 
frared images indicate a central cavity. Instead this emission must 
arise for a component which reaches higher temperatures, either 
because it consists of smaller grains, or receiving additional heat- 
ing. This component has a much smaller mass than the cold dust. 

The extinction allows us to map the cold dust at much higher 
resolution than would otherwise be possible. Imaging thermal 
emission of cold dust requires far-infrared space telescopes such as 
ISO and Spitzer; their small apertures are "resolution-challenged" 
(Rie ke et al.l20oJ) : e.g. the 6-arcsec resolution at 24/im for Spitzer 
would leave the core of NGC 6537 unresolved. In contrast, our ex- 
tinction map has a resolution of 0.14 arcsec per pixel. As the dust 
shell closest to the core was ejected last, this gives a unique han- 
dle on the morphology of the final AGB wind, and helps constrain 
the origin of the bipolar morphology. 

The extinction map is far more symmetric than either the 
HST or the infrared images. The structure delineated by the in- 
ner cavity shows an elongation approximately (but not perfectly) 
along the bipolar axis, with an axial ratio of 1.5:1. For compar- 
ison, the optical image shows an elongation of the full structure 
(at the faintest emission levels) of approximately 4:1. A large am- 
plification has clearly taken place. The structure shows that the 
mass loss on the AGB was much more spherically symmetric than 
the current PN would suggest. 

The shell is relatively thick, with an outer radius at least 
twice the inner radius. The optical emission appears embedded 
in the dust shell, and so the extinction map measures the extinc- 
tion from the front shell only (this is the reason we can apply the 
extinction map for the central star). The observed ratio of peak 
extinction over central extinction (in mag) of 1.4, would require 
an outer radius of the dust shell of 6 arcsec, for an observed inner 
radius of 3 arcsec and uniform density. However, the extinction 
includes an ISM foreground component, and the observed extinc- 
tion ratio is therefore a lower limit. A higher ratio corresponds to 
a smaller outer radius. 

The dust shell may well extend further out than derived 
from the extinction map: the optical H/3 flux becomes too faint 
to accurately determine extinction values. But the estimate above 
makes is likely that the dense part of the shell does not extend 
much beyond the limits of the map. The AGB wind will extend 
much further asal/r 2 envelope, as suggested by the large extend 
of the optical emission. But an extended 1/r 2 envelope would not 
give rise to a large difference in extinction over the (small) central 
hole. We conclude that inner dust shell is distinct from the outer 
AGB wind, either tracing a short lived (< 10 4 yr) event of very 
high mass loss, and/or swept-up gas. 

The final point to note is that we c an exclude the 'thin' 
almost-edge-on disk which was proposed bv lCuesta et al.l f.995). 
This would give rise to a narrow, extinction structure, with a 
strongly elongated 'higher extinction region' perpendicular to the 
polar axis. The observed structure is only mildly elongated, and 
the direction of elongation is opposite to what would be expected 
from a thin disk. The observed elongation is reminiscent of a 
cylindrical structure, or of a torus with thickness comparable to 
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Figure 13. Schematic view of shell. Circumstellar shell may be 
elongated along the holes so as it looks like a 'torus'. 

its radius (Tig. 1131 . but apart from the polar holes the structure 
show no clear deviation from spherical symmetry. 

An inevitable question is whether this dust shell is respon- 
sible for collimating the bipolar flow. In favour of this is the lo- 
cation, at the neck of the flow, and the compact size, 10 per cent 
of that of the lobes. The similarity of dynamical age between the 
dust shell and the shell shows that the collimation first occurred 
very early in the evolution, during the AGB or in the early post- 
AGB wind. This may argue against versions of the Interacting 
Stellar Winds where the shaping is due to the fast wind from the 
central star of the PN, as this wind only develops in the later 
post-AGB evolution, at higher stellar temperature. 

The timing of the lobes is consistent with current ideas that 
the shaping originally occurs due to a change in t he wind mor- 
pholo gy at the end of the high mass-loss phase (jZiilstra et al.l 
2001). A model which could fit the observed characteristics (a 
symmetric she ll with lower d ensities only towards the poles) is 
presented bv lSahai and TraugeJ (l998): they propose that a fast 
jet-like wind punches holes in the spherical AGB shell, and that 
the subsequent post-AGB evolution amplifies this sudden depar- 
ture from spherical shape. 

Such a model may require a binary companion as source 
for the jet. There is no evidence either against or in favour of 
a binary at the core of NGC 6537. The optical magnitude rules 
out a solar-like star, but a fainter red dwarf or a white dwarf 
are not excluded. The near infrared detection may show some 
excess emission, but this cannot be taken as evidence for a red 
companion. For comparison, the central star of the very similar 
PN NGC 6302 has been show n to have an infrar ed excess which 
is not a stellar continuum (Matsuu ra et alj l2005V 



4.3 Hierarchical structure 

In NGC 6537, a hierarchical structure is found, with the peak 
of PAH intensity located further out than that of Ha. This 
is consistent with obs ervations in oth er objects (NGC 70 27 
dWoodward et al.lll989ft and NGC 6302 dMatsuura et al.ll2005ft N ). 
where PAH emission is located outside of the ionised gas. The 
extinction continues to increase outside the PAH region, and 
PAH emission is not found where the extinction is high (~5 arc- 
sec). This can be a consequence of the lack of photo excita- 



tion in the high exti nction region, or a difference in abundances 
dMatsuura et al. 112004ft . 

The abundance analysis from optical line ra tios shows that 
NGC 653 7 is oxygen-rich {Pottasch et al.ll2000bft . The presence 
of PAHs (Roche and Aitken 1986), which are normally assumed 
to be formed from carbon-rich gas, is unexpected but this chem- 
ical d ichotomy is also seen in a few other PNe ( Mat suura et all 
2005). PAHs have smaller absorption coefficients in the optical 
than silicate grains, and PAHs are unlikely to be the source of 
high extinction in the optical: amorphous silicate is likely the 
main contributor to the optical extinct ion. Crystalline sili cates 
have also been detected in NGC 6537 (|Molster et alJ 2002). As 
the peak of extinction occurs outside of the PAH emission peak, 
the silicate grains are located outside the location of the PAHs. 
This can indicate that the silicates pre-date the PAHs. 

Hydrogen molecules are detected in NGC 6537: the brightest 
H2 region is found a round the ionized core, corresponding to the 
high extinction shell (jDavis et al.l2003^ . The presence of H2 shows 
that the dust shell is partly molecular, and the separation between 
extinction and Ha traces an ionization front. The PAH emission 
is located close to the ionization front but the precise relation, 
and whether the PAHs formed in-situ or also date from the AGB, 
remains to be determined. 

The hierarchy, with the PAH located in between the ionized 
and molecular gas, suggests that the PAH emission may be as- 
sociated with the photo-dissociation region of the molecular dust 
shell. Following (jMatsuura et al.l2004ft . we associate the hot dust 
component (Section 14.21 with dust located in the PAH-emitting 
region. 



5 CONCLUSION 

The high-resolution images presented in this paper detect the cen- 
tral star of NGC 6537. The agreement between the optical and 
infrared images secure the identification. The magnitudes are rel- 
atively uncertain, due to the high background, but the derived 
values agree within the uncertainties with a hot black body. The 
luminosity depends on the assumed temperature: dynamical con- 
siderations, comparing with the stellar evolution tracks, argues 
for relatively low luminosities ~ 1O 3 L0. This gives a temper- 
ature in the range 1.5-2.5xl0 5 K, slightly lower than the Hell 
Zanstra temperature of 3.4 x 10 5 K. We derive a core mass around 
M c w 0.7-0.9 M , and a progenitor mass Mi = 3-7M . 

The progenitor mass is consistent with the classification of 
NGC 6537 as a type-I PN. Bipolar PNe have been suggested to 
have higher than average initial mass: our results do not conflict 
with this. 

The extinction map is not uniform inside NGC 6537, but 
shows a shell structure of inner radius of 2arcsec and outer ra- 
dius of at least 4arcsec. The central star is located at the cen- 
tre of the near-symmetric extinction distribution. The dust shell 
traces the high mass-loss phase during the AGB phase, and was 
expelled from the central star about 1400 years ago. There is no 
evidence for a thin disk, as has been suggested. The dust shell 
is far more compact, and far more symmetric than the extended 
bipolar lobes. A thick torus or cylindrical structure can fit the 
observed extinction map. 

The dynamical age of the lobes and dust shell are very simi- 
lar. This suggests that the origin of the bipolarity should be close 
to the end of the AGB. The largely symmetric dust shell, com- 
bined with the strongly bipolar lobes, are consistent with models 
of a fast change of wind structure at the end of the AGB mass 
loss. 

The structure shows radial hierarchy, with the ionized core 
located at the inner edge of the dust shell, and the PAH emission 
in between, tentatively associated with the photo-dissociation re- 
gion. 
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